DETERMINANTS OF WEIGHTED PATH MATRICES 



KELLI TALASKA 

Abstract. We find rational expressions for all minors of the weighted path matrix of a di- 
rected graph, generalizing the classical Lindstrom/Gessel-Viennot result for acyclic directed 
graphs. The formulas are given in terms of certain flows in the graph. 



1. Introduction 

One of the most elegant results in algebraic combinatorics is that expressing minors of the 
weighted path matrix of an acyclic graph as a sum over all collections of pairwise vertex- 
disjoint paths connecting the specified sources and sinks. While this result was (re) discovered 
and popularized in combinatorics by Gessel and Viennot j GV85j . the key ideas can also be 
found in an earlier, more abstract, formulation by Lindstrom jLin73j and in Karlin and 
McGregor's work on coincidence probabilities |KM59] . 

In the acyclic setting, the weighted path matrix has entries which are polynomial in the 
edge weights of the graph. The main result of this paper is an extension of the above result 
to graphs which are not necessarily acyclic; in this setting, the weighted path matrix has 
entries which are formal power series in the edge weights. In particular, for each minor, we 
give a rational expression whose numerator and denominator are both polynomials given by 
simple combinatorial formulas. 

The key result in this paper is closely related to work in several other papers. First, we 
note that Viennot 's work on the theory of heaps of pieces |Vie89j can be used to derive 
our rational expressions for the entries of the weighted path matrix, but it does not appear 
that this approach can be used to obtain the formulas for larger minors. Fomin's work on 
loop-erased walks |Fom01] also examines the weighted path matrix, giving an expression for 
each minor as a sum indexed by an infinite but minimal collection of path families satisfying 
certain intersection criteria. The path families indexing Fomin's formulas depend on a choice 
of labeling of the sources and sinks, but our indexing path families are uniquely determined. 

The proof of our result involves an involution which generalizes the "tail-swapping" proof 
from the acyclic case. The involution is adapted from the author's work on total positivity in 
Grassmannians |Tal08j . in which planar graphs with directed cycles play a key role, though 
the analogues of weighted path matrices are rather different in the Grassmannian setting. 
Besides planarity, there are a number of topological and other technical conditions placed 
on the graphs used in |Tal08] . none of which are necessary in this paper; our formulas hold 
for all directed graphs. 

The classical result for acyclic graphs has many applications, including combinatorial 
proofs of the Jacobi-Trudi determinantal formulas for Schur functions and MacMahon's 
formula for the number of plane partitions (see |GV85j and [AigOl] for these and several more 
well-known examples). We hope that, with the added flexibility of working with directed 
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cycles, our generalization will have many applications as well. One such example comes from 
algebraic statistics - this paper was partly motivated by the author's work with Draisma 
and Sullivant on Gaussian graphical models. A forthcoming sequel to [STDlQj will use the 
results of this paper to address models with directed cycles. 

2. Statement of the main theorem 

Let G be a directed graph. Loops and multiple edges are permitted. We impose no further 
conditions on G. In particular, G is permitted to have directed cycles, and we do not assume 
that G is planar. Let V = {vi, . . . , Vn} be the vertex set of G. Assign to each edge e of G 
the formal variable Xe', we call Xe the weight of the edge e and assume that all edge weights 
commute with each other. 

A path P = (ei, 62, . . . , em) in G is formed by traversing the edges Ci, 62, . . . , e^, in the 
specified order. We write P : v --^ v' to indicate that P is a path starting at a vertex v and 
ending at a vertex v'. 

Define the weight of a path P = (ci, . . . , em) to be 

Wt(P) = Xe,---Xe^. 

A path P : V V with no edges is called a trivial path and has weight 1. Paths may or 
may not have self-intersections; those with no self-intersections will be called self- avoiding. 

Definition 2.1. The weighted path matrix of G is the matrix M whose entries mij are the 
formal power series 

"^ij = ^ wt(P), 

the sum over all directed paths P : Vi vj. 

This generalizes the classical definition for acyclic graphs. Note that if A is the weighted 
adjacency matrix of G, then M = (J — A)~^. 




V2* 



V3 



Figure 1. An example of a graph G with a directed cycle. 



Example 2.2. Consider the submatrix Mji 2},{3,4} of the weighted path matrix of the graph 
G in Figure [T| For the matrix entry rrii^, we note that any directed path from Vi to ^3 must 
traverse the path with weight abdfh, and it may complete an arbitrary number of loops of 
the cycle of weight bde along the way. Similarly, for mi4, any directed path from vi to ^4 
must traverse the path with weight abc, and it may complete an arbitrary number of loops 
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of the cycle of weight bde. There is a unique path from V2 to V3, with weight gh, and there 
are no paths from V2 to V4. Thus, we have 



Since this particular example is fairly simple, we can easily see how to write the entries of 
the matrix as rational expressions: 



The objective of this paper is to provide such rational expressions for the entries and larger 
minors of the weighted path matrix when the network may be more complex. We now define 
several families of paths and cycles needed for the statement and proof of the main theorem. 

Definition 2.3. Let C{G) denote the set of collections C of self-avoiding cycles in G which 
are pairwise vertex-disjoint. Each cycle in a nonempty collection C much contain at least 
one edge, and we consider the empty collection an clement of C{G), with weight 1. 

For the remaining families, let A = {ai, . . . , a^} and B = {bi, . . . , bk} be two subsets of 
the vertex set V, and let tt be a permutation in the symmetric group Sk- 

Let Va,b,tt{G) denote the set of path collections P = (Pi, . . . , Pk) such that Pj is a directed 
path in G from vertex Oj to vertex b^^i^i). Set Va,b{G) = 1) Va,b,-w{G). 




Let Va,b{G) be the subset of path collections in Va,b{G) such that 

• each Pi is self-avoiding, and 

• the paths Pj and P, are vertex-disjoint whenever i ^ j. 

Let J-a,b{G) denote the set of self-avoiding flows connecting A to B, i.e., pairs (P, C) 
such that 



• P and C are vertex-disjoint. 

Informally, a self-avoiding flow in is a collection of self-avoiding paths connecting 

A to B along with a (possibly empty) collection of self-avoiding cycles such that the paths 
and cycles are all pairwise vertex-disjoint. 

Definition 2.4. Define the weight of any path collection, cycle collection, or self-avoiding 
flow to be the product of the weights of the paths and cycles it contains, counting the weights 
of any repeated edges with the appropriate multiplicities. 

Define the sign of a path collection, cycle collection, or self-avoiding flow as follows. If 
P £ 'Pa.b,tt{G), set sgn(P) = sgn(7r). If C G C{G), let |C| denote the number of cycles in 
C, and set sgn(C) = (—1)'^'. Each self-avoiding flow F has a unique decomposition as the 
disjoint union of a path collection P and a cycle collection C; set sgn(F) = sgn(P) sgn(C). 






. P G VaAG), 
• C G C{G), and 
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Theorem 2.5. Suppose G is a directed graph with weighted path matrix M . Then the minor 
Aa,b{M), with rows indexed by A and columns indexed by B, is given by 

sgn(F) wt(F) 

(2.1) A,MM) = 

sgn(C)wt(C) 

ceC(G) 



In general, the fraction in equation (2.1) is not reduced, i.e. the numerator and denom- 
inator share a common factor. Consider the connected components Hi,H2, . . . ,Ha of the 
subgraph Go consisting of all edges which appear in at least one directed cycle of G. It can be 
easily verified that the denominator above can be factored as a product over the connected 
components of Go, i.e. 

a 

sgn(C)wt(C) = n E sgn(C)wt(C). 

CeC(G) i=l C£C(Hi) 

Each factor sgn(C)wt(C) corresponding to a connected component Hi of Go is irre- 

ducible, and it will cancel if and only if every collection of (self-avoiding, pairwise vertex- 
disjoint) paths in Va,b{G) avoids the component H^. 



Example 2.6. Let us use Theorem 2.5 to compute some minors of the weighted path matrix 
of the graph G in Figure [l} There are two cycle collections in C(G), the empty collection and 
the single cycle of weight bde, so the denominator of every minor, before canceling common 
factors, will be 1 — bde. 

The matrix entry mi3 is the 1x1 minor A|i}^|3}(Af ). There is a unique self-avoiding flow 
from vi to f3, namely the path of weight abdfh. Thus, equation (2.1) tells us that 

abdfh 



1 — bde ' 



which is consistent with our calculations in Example 2.2 



For the matrix entry = A|2},{3}(M), we again have a unique self-avoiding path, but 
in this case there are two self-avoiding flows, since the cycle of weight bde does not intersect 
the path of weight gh. Thus, 

gh + gh{—bde) gh{l — bde) ^ 
1 — bde 1 — bde ' 

as we previously showed. 

If we look at the 2x2 determinant |M|i 2},{3,4}|, we find a single self-avoiding flow connect- 
ing A = {ai = Vi,a2 = V2} to B = {bi = V3, 62 = ^^4}, and the sign of this flow is negative, 
since the pair of paths connecting AtoB sends ai to 62 and 02 to 61. Equation (2.1 ) implies 
that, as expected, the determinant is 

—abcgh 



|^{l,2},{3,4}| 



1 — bde 



Proof of Theorem 2^: Using the Leibniz expansion of the determinant, we have 
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Using Definitions 2.1, 2.3, and 2.4, we can rewrite this as 



Aa,b(M) 




sgn(7r) wt(P) = J2 sgn(P)wt(P). 

Thus, to prove Theorem |2.5[ it would suffice to show that 

^ J2 sgn(P)wt(P)sgn(C)wt(C)= ^ sgn(F)wt(F). 

ceC(G)pep^_5(G) FeJ-A,s(G) 

That is, all terms on the left cancel except those for which (P,C) G J^a,b{G), since 
each self-avoiding flow F = (P, C) satisfies sgn(F)wt(F) = sgn(P) wt(P) sgn(C) wt(C). 
We prove this by constructing a weight-preserving and sign-reversing involution on pairs 
(P,C) ^J^aAG). 

Suppose that P = (Pi,...,Pfe) G VaM^) and C G C{G), but (P, C) ^ J^aA^)- This 
means at least one of the following must be true: 

• there exists an i such that Pi is not self-avoiding, 

• there exists a pair i i' such that Pi and Pj' share a common vertex, or 

• there exists an i such that Pj and C share a common vertex. 

Define (f(P,C) = (P*,C*) as follows. Choose the smallest i such that Pj is not self- 
avoiding or shares a common vertex with C or with some Pj/ with i' > i. Then, following the 
algorithm below, we either swap the tails of two paths, as in Figure |2} or we move a cycle 
from Pj to C or vice versa, as in Figure |3| 

Let Pj = (ei, . . . , em)- Let tail(ej) denote the tail of the edge Cj. (Each edge points from 
its tail to its head.) Choose the smallest q such that the vertex tail(eq) lies in C or in some 
Pi' with i' > i, or tail(eg) = tail(er) for some r > q. 

• If tail(eq) lies in some Pj/ with i' > i, choose the smallest such i'. (This case allows 
for the possibility that Pj also intersects itself or C at the vertex tail(eg).) We then 
swap the tails of Pj and Pj/ as follows. Let Pj/ = {hi, . . . , hm'), and choose the 



smallest q' such that tail(/ig/) = tail(eq). Set Pj* = (ei. 



-ihq'hn 



hr 



and P*, = {hi, /iq/_i, eg, e^+i, . . . , e^). Set P* = P \ {P„ Pj/} U {P* , P*,] and set 
C* = C. Note that P* ^ P, since Pj and Pj/ have different endpoints. Q 




i > 




Figure 2. Tail swapping. The paths Pj and P* are drawn as solid curves, 
and the paths Pj/ and P*i are drawn as dashed curves. 



Technically, we have swapped the heads of the two paths, but it is common to refer to such an operation 
as tail swapping. 
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• Otherwise, we find the first point along Pi where we can move a cycle from C to 
Pi or vice versa, as follows. If Pi is not self-avoiding, let ^ be the first cycle that 
Pi completes. More precisely, choose the smallest t such that tail(es) = tail(et) for 
some s < t; then i = {es,es+i, ■ ■ ■ ,et_i). If Pi is self-avoiding, then set t = oo. If 
C intersects Pi, choose the smallest u such that tail(eu) appears in a (necessarily 
unique) cycle L = {li,l2, ■ ■ ■ Jw) in C, where tail(/i) = tail(e„). If C and Pi are 
vertex-disjoint, then set u = oo. At least one of t or u must be finite, and t ^ u. 
Note that both t and u may be distinct from g, as is the case in Figure [3j 

o lit <u, move i from Pi to C. Set C* = C U P* = (ei, . . . , e^-i, e*, . . . , e^), 

and P* = P \ {Pi} U {P*}. 
o If u < t, move L from C to Pi, in the following position. Set C* = C \ {L}, 
P* = {ei,..., e,_i, h,..., L, e,, . . . , e^), and P* = P \ {P,} U {P*}. 




Figure 3. Moving a cycle. The paths Pi and P* are drawn as solid curves, 
and the cycle collections C and C* are drawn as dashed cycles. 

It is easy to see that, with this definition, the image (P*, C*) is again a pair of the required 
kind, i.e. P* e Va,b{G), C* G C{G), and (P*,C*) ^ J^a,b{G). 

Let us verify that ip is an involution. First, we check that if does not change the value of 
the index i. That is, among all paths in P* which have self-intersections or intersections with 
other paths in P* or cycles in C*, the path with the smallest index is P*. Indeed, our moves 
only affect Pi, Pi', and C, keeping their combined set of edges intact, so the involution will 
not introduce a new self-intersection in any path Pj such that j < i, nor will it introduce a 
new intersection between such a path Pj and any other path or cycle. 

Consider the first case, in which ip(P, C) = (P*, C). After swapping tails, P* still has no 
intersections with C or any of the other paths before the vertex tail(eg). Further, P* does not 
have any self intersections before tail(eg) - though it may have self-intersections at tail(eg) - 
since Pi did not have any self-intersections before tail(eg) and the tail of Pi' did not intersect 
Pi before tail(eg). Thus tail(eq) remains the first vertex along P* with an intersection. Now, 
P* and P*, intersect at this vertex, and no path with smaller index intersects P* at tail(eg), 
so applying ip again swaps the same tails. 

Consider the second case, in which !f{P, C) = (P*, C \ {L}) or if(P, C) = (P*, C U {£}). 
Here Pi intersects itself or C at the vertex tail(eq), but does not intersect any other path 
at this vertex. After moving a cycle, the same is true for P*. If the cycle moved starts at 
tail(eg), i.e. tail(eg) = tail(ei) or tail(eg) = tail(e„), then either a self-intersection becomes 
an intersection with C, or an intersection with C becomes a self-intersection. If the cycle 
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moved starts later, then the intersections at tail(eg) remain the same type after applying if. 
If Pi intersects C before completing its first cycle, then P* will complete its first cycle before 
intersecting C \ {L}. If Pj completes its first cycles ^ before intersecting C, then P* will 
intersect C U {£} before completing its first cycle. Thus, the same cycle is moved twice if we 
apply </? twice, which completes the proof that tp is an involution. 

It remains to verify that ip is sign-reversing and weight-preserving. In the first case, 
with (^(P, C) = (P*,C), the permutations corresponding to P and P* differ by a single 
transposition, so sgn(P*) = — sgn(P). The path collections P and P* use the same multiset 
of edges, so wt(P*) = wt(P). Since C remains fixed, sgn(C) and wt(C) introduce no changes. 

In the second case, with (^(P, C) = (P*, C\{L}) or (^(P, C) = (P*, CU{£}), the endpoints 
of paths in P are the same as those of the corresponding paths in P*, so sgn(P*) = sgn(P). 
However |C*| = |C| ± 1, so sgn(C*) = — sgn(C). Though edges are moved between P 
and C, the pairs (P, C) and (P*, C*) both use precisely the same multiset of edges, so 
wt(P*) wt(C*) = wt(P) wt(C). 

In each case, the overall sign changes, and the weight is preserved, completing our proof. 

□ 
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